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1. Summary 

This report conti ins the results of the first noise susceptibility computer 
simulation tests of the complete MIA receiver analytical model. Computer 
simulation tests have been conducted with both Gaussian and pulse noise inputs. 

The results of the Gaussian noise tests have been compared to results predicted 
previously (Ref. JSC-09083) and have been found to be in substantial agreement. 

The results of the pulse noise tests will be compared to the results of planned 
analogous tests in the Data Bus Evaluation Laboratory at a later time. The MIA 

computer model is considered to be fully operational at this time although 
refinements and modifications will be made as further information on the 
Singer MIA design becomes available. 

2. Introduction 

This report has been prepared to present the initial results from the 
KLA analytical model and to act as an interim status report. 

3. Discussion 

A. The computer model of the Singer MIA receiver has been described 
in a previous paper (Ref. 1.3-DN-C0104-002). In that paper the MIA 
filter had been described in general terms, but had not been 

implemented in the computer model. The entire MIA simulation program, including 
the MIA filter has now been implemented and is described below. 

B. Decoder model and input circuitry 

1. The decoder section of the computer model has not been changed. 

(This section is described in Ref. 1.3-DN-C0104-002). 

2. The clipper level has been chosen to be + 16 volts. This 
value has been chosen to minimize the effects of high amplitude 
noise while allowing the largest permissible signal voltage 

(+ 15 volts) to pass unattenuated. The actual clipper level used by the 
Singer design is still unknown. 
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3 . The transformer and threshold detector are as described previously. 
(Kef. 1.3-DK-C0104-002i 
C. Filter Models 

1. Two filter models have been developed using bi-linear Z 
transform techniques. (See Appendix A). A 6-pole Butterworth low-pass 
filter and a 6-pole Bessel low pass filter have been implemented. 

2. The Butterworth filter has been tested using sine 

wave and step response inputs. The results are shown in Figures 1 
and 2 respectively. These results indicate that the filter is 

• functioning as desired. 

3. Th9 Bessel filter has also been tested using sine wave and 
step response inputs. The results are shown in figures 3 and 4. 

These results also agree with expected results . 

4. At the time of the implementation of the MIA simulation program, 
the type of filter being used by the MIA was uncertain. The first 
filter modeled was of the Butterworth type. The Bessel filter was 
modeled in order to make comparisons with the results of a previous 
analysis of data bus noise effects (Ref. JSC-09083). The Gaussian 
tests were run using the Bessel filter, while the pulse tests were 
run using both types of filter. 


4. Results 

A. Gaussian Noise tests 

1. The Gaussian noise tests were undertaken as an attempt to 
verify the operation of the MIA cor niter model. A paper by George 
Proch of LEG (Ref. JSC-09083, referred to hereafter as the JSC paper) 
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was prepared which predicted the performance of the MIA with Gaussian 
noise input. It was desired to run simulations with the MIA model and 
compare the results with the predictions in order to confirm the accuracy 
of both the JSC paper and the MIA model. 

2. The results of the Gaussian tests using the Bessel filter 
are plotted in figure 5 together with the predictions from the 
JSC paper. The probability of an incorrect word transmission 
vs RMS noise amplitude is the quantity compared. Other * 
quantities predicted by the JSC paper (such as undetected word errors 
or bit decision errors) result in numbers which would require 

an inordinate number of computer runs and vast amounts of 
computer time. Two sets of curves are plotted for the JSC 
results. The dashed lines plot the original predictions. 

It was subsequently determined that the equivalent noise 

bandwidth of the filter analysed by the paper should be approximately 

1.06 x F c and not 1.66 x F c as calculated originally. This has been 
confirmed by Mr. Proch, and the solid lines represent the revised data. 

3. As shown in figure 5, the results from the MIA model tests 

agree quite closely with the dependent sample case from the JSC paper. 

The following information should be taken into consideration when 
comparing the results: 

a. The JSC prediction should be considered to be optimistic. 

The assumptions used to get the results (square input data 
plus Gaussian noise) have been undercut by the recalculation 
of the step-response rise time of the input filter. The 
original assumption made was that the rise time was fast 
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enough to allow the assumption of square data input. A 
comparison of figure 5 in this report with figure 6 in the JSC 
report indicates that the rise time of the filter is over 
2 times slower than predicted by the JSC paper. (This slower rise 

time has been confirmed by Mr. Proch in a lab test of an analog 
Bessel filter.) This new rise time will tend to round 
the input data and cause more "no decision" samples than 
predicted . 

b. The JSC paper predicted that the actual results would fall somewhere 
between the dependent and independent sample cases which were 

to be considered as limiting cases. It is felt that due to 
the fairly slow rise time of the filter, it would be 
expected that the actual results would fall closer to the 
dependent sample case since the low-pass filter acts to 
limit the change between samples. 

c. Gaussian noise was simulated by using a random number generator 
which produces a Gaussian distribution with zero mean and a 
variance of one. These numbers are generated at a simulated 

rate of 96 megasample3 per second. This gives a band-limited noise 
spectrum which is white to 48 megahertz. In order to conform 
to the JSC paper's assumption of a 4 megahertz noise bandwith, 

the RMS value of the generated noise was multiplied by 
or .289. 

4. It is felt that the comparison of the JSC paper predictions and the 
MIA model test results, together with the detailed analysis which 
went into the resolving of initial discrepancies, has resulted in a 
high level of confidence in both the JSC paper (as modified) and the 
MIA computer model. 
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B. Pulse tests 

1. A series of pulse tests have been run with the MIA computer 
model. Each test consists of the reception of a specified 
number of 28 bit words with one noise pulse added per word.. 

The noise pulse is positive polarity, and its position in the 

word, is selected at random. Each subseries of tests consists , 

of adding pulses of constant width while varying the pulse 

amplitudes. The results obtained vising the Bessel filter are 

plotted in figure 6, and those obtained using the Butte rworth 

filter are plotted in figure 7. In both cases the quantity 

plotted is the probability (in percent) of a word error vs 

i 

the pulse amplitude in decibels relative to one volt. Pulse 
widths are integer multiples of 10.4 nanoseconds (resulting from 
the 96 megasample/second rate). 

2. The results in both tests are qualitatively as expected. For 
very short pulse widths, the pulse amplitude required to cause 
a word error is high due to the filter attenuating the pulse. 

As the pulse width increases, the amplitude needed to cause a word 
error decreases. For very long pulse widths, the probability of a word 
error approaches 100 % for any amplitude above the error threshold 
level (approximately 1 volt). 

3. The results of the Bessel and Butterworth tests are plotted 
together in figure 8 to compare the results in the two cases. 

It is apparent that both filters operate similarly at high 
noise levels. However, for pulse amplitudes close to the threshold 
level it can be seen that the Bessel filter performs better in all 
cases. This is attributable to the better phase characteristics 
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of the Bessel filter. A comparison of the step responses of the 
two filters (figures 2 and U) shows that the Butterworth filter 
exhibits considerable overshoot and ringing resulting in substantial 

degradation of the data wave form. The Bessel filter 
exhibits virtually no overshoot and the waveform is less 
distorted. 

3. The results of these pulse tests will be compared with 

analogous tests to be carried out in the Data Bus Evaluation 
Laboratory. These two sets of data will then be used in 

conjunction with a data bus impulse noise model currently being 
developed to determine error probabilities on the data bus lines. 

5 . Conclusion 

This report has presented the results of the initial computer simulation runs 
using the complete MIA receiver model. Results have been given for both 
Gaussian and pulse type noise and for both Bessel and Butterwcrth type input 
filters. Based on these results, the MIA computer model is considered fully 
operational at this time. The results illustrated in Figures 5-b provide a 
basis for the evaluation of the. results of planned tests on the MIA in the 
Data Bus Evaluation Laboratory. It appears that a Bessel type input filter has a 

slight advantage over a Butterworth filter for pulse type noise. It should be 
noted that in no case did a noise pulse of peak amplitude less than 1 volt 

cause a word error. This indicates that only noise waveforms with peak 
amplitudes greater than or equal to 1 volt will be capable of causing 
errors on the data bus lines. 
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The results presented in this report are tentative since the MIA computer 
model is based on preliminary information (albeit the most current available). 
Further updates to the computer model will be made as information becomes 


available 





6-pole Eutterworth Filter Step Response (f c 1.5 MHz) 

Figure 2 
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1, Bessel Filter 

The transfer function for a 6-pole Bessel (naximally flat delay) filter 
is veil-known to be: 


°12 (S> 


10.395 


10,395 + 10,395Ts ♦ l* t 725?2s 2 + l,260? 3 s 3 + 210T 1 *s l * + 21T 5 s 5 ♦ T 6 s 6 
where T is the desired delay in seconds. Using Fig, 13-16, pg. 391 
of Modern Network Synthesis (Ref. M. E. Van Valkenburg, Modern Network 
Synthesis . New York: Wiley & Sons, 1967), for a 6-pole Bessel filter 

the 3-db down point corresponds to 
cj c T = 2.7 
T = 2 *7 


Using the bi-linear Z transform requires that we set 

w c . DT 
coa = tan — 2 


to compensate for a frequency axis warping introduced by the transform, 
where DT is the sample interval, and acts as a dummy variable. 

We now set 


W A 

T ’-srr 


A = 

10,395 

B = 

10.395T 

C = 

14.725T 2 

D = 

i,26ot 3 

E = 

210T 1 * 

F = 

21T 5 

G = 

T^ 


Letting, 
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imd performing the substitution 
2-1 


S « 


2 + 1 


which is the bi-linear Z transformation, we obtain 

, v _ A 

G 12 (») 




r- - 1? + Ei 

( z '~ A 


r* - 

} ♦ G/2-S 

A + B (z ♦ l) + C| 

(tvtJ * 1 

ITTT? ^ 

lz * l) 

lz ♦ 1^ 

'2 ♦ 


Multiplying through by (2 + 1)^, cna letting 
Q1-A + 3 + C + D + E + F 
Q2 * 6A ♦ UB = 2C - 2E - IF- 6G 

Q3 » 15A ♦ 53 - C - 3D - E + 5F ♦ 15G 

QU - 20A - UC ♦ 4E - 20G 

Q5 ■ 15A - 53 - C - 3D - E - 5F ♦ 15G 

Q6 » 6A - 4B + 2C - 22 ♦ 4F - 6G 
Q7-A-3+C-C+E-F+G 
We obtain 

/v A (Z 6 + 6 Z 5 ♦ 15z U + 20t 3 ♦ 15z 2 + 62 + 1) 

G 12 ' Z ‘ = — i 1 3 2 

Qlz^ ♦ Q 2 z' + Q 3 z* ♦ Q^z + Q 5 z + Q 6 z ♦ Q7 

b A (1 ♦ C7T 1 * 15z" 2 ♦ 20 z “ 3 + 13Z" 14 + 62 5 ♦ 2~ 6 ) 

Q1 ♦ QBz" 1 + ^3 -2 ♦ ^z-3 + Q5z*"^ + Q6 z“^ + Q7z ^ 

Since the z " 1 operator corresponds to a delay of one sample tine we finally 
obtain 

Y (nT) = -3- [ A [ X(nT) + 6 X(nT - T) - 15 X(nT - 2T) ♦ 20 X(nT - 3T) 

Q1 

+ 15 X(n? - IT) ♦ 6 X(nT - 5T) + X(nT - 6 t) ] 

- (Q2 Y(nT - T) ♦ Q3 Y(nT - 2T) + QUY(nT - 3T) 

♦ Q 5 Y(nT - 4T) ♦ Q 6 Y(nT - 5T) ♦ Q7 Y(nT - 6 T)]] 

Where X(nT) is the sampled input signal and Y(nT) is the sampled output 
(filtered) signal. 
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2. Buttcrvorth Filter 

The transfer function for a 6-pole Buttervorth (maximally flat amplitude 
response) filter is: 


°12 <»> 


w c 6 


y p 6 ♦ 3.6637 s w „^ ♦ 7.^641 s^u ** ♦ 9.ll*l6s\> ^ ♦ 7.1*61»ls 1 * w ♦ 3.8637s w 

» c c c » c 


Where u is the 3-db cutoff frequency, 
c 

Using the bi-linear Z transform requires (see Bessel Filter, Part 1) 
that ve set 


! 


w. 


U> • DT 
tan ■ . c ■ , 


LettinG 


Ql » 3.8637 w A 
Q2 * 7.1*61*1 u A ^ 
Q3 * 9.1^16 a* 3 
Ql* * 7.1*61*1 uj^l* 

Q5 * 3.8637 w A 5 
Q6 = uj* 


We have 
G 12 ( 5 ) ■ 


06 


s 6 ♦ Qls 5 ♦ Q2s U + 03s 3 + Ql*s^ ♦ Q5s + Q6 


Taking the bi-linear Z transform, multipying by (z ♦ 1) D and letting 

Cl = 1 ♦ Ql ♦ Q2 ♦ Q3 ♦ QU ♦ Q5 ♦ Q6 

C2 = - 6 - 1*Q1 - 2Q2 ♦ 2Ql* ♦ 1*Q5 + 6 q6 

C3 * 15 ♦ 5Q1 - Q2 - 303 - Ql* ♦ 505 ♦ 1506 

CU = - 20 ♦ 1*Q2 - l*Ql* ♦ 20Q6 

C5 “ 15 - 501 - 02 + 2Ql* - 1*Q5 ♦ 606 

C6 ■ - 6 ♦ 1*01 - 2Q2 ♦ 2QU - 1*05 ♦ 6q6 

C7 = 1 - 01 ♦ 02 - 03 ♦ Ql* - 05 ♦ Q6 


MV 
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„ , v C6 (z 6 ♦ 6z 5 ♦ I5z u ♦ 20s 3 ♦ 15z 2 ♦ 6z ♦ 1) 

c 12 (z) * 

Clz 6 ♦ C2z 5 ♦ C3Z 1 * ♦ CUz 3 + C5z 2 + c( ; z + c7 

, C6 (It- Cz" 1 » '5z~ 2 ♦ 20z~ 3 + 15z~ U + 6z" 5 ♦ z~ 6 ) 

Cl ♦ CCz" 1 ♦ C3z“ 2 ♦ CUz” 3 ♦ C5z”^ ♦ C6z" 5 ♦ CTz“^ 
which yields 

Y(nT) «-JL [ [Q6 ( X(nT) + 6 X(nx - T) ♦ 15 x (nT - 2T) + 20 X( n T -3T) 

Cl 

♦ 15 X(nT - UT) ♦ 6 X(nT - ST) ♦ X(nT -6T) 1 

- C2 Y(nT - T) ♦ C3 Y(nT - 2T) ♦ C4 Y(nT - 3T) 

♦ C5 Y(nT - 1*T) ♦ C6 Y(nT - 5T) ♦ C7 Y(nT - 6T)]] 

Note 

The implementation of each of the above equations is of the "direct 
implementation" fora. This is the easiest fora to implement, but it does have 
some disadvantages. The primary problem is that the coefficients, necessarily 
being of finite length, begin to cause erroneous output when w c ^ T becones too 


small. 



